Monodisperse domains by proteolytic control of the coarsening instability 
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The coarsening instability typically disrupts steady-state cluster-size distributions. We show that 
degradation coupled to the cluster size, such as arising from biological proteolysis, leads to a novel 
fixed-point cluster size. Stochastic evaporative and condensative fluxes determine the width of the 
fixed-point size distribution. At the fixed-point, we show how the peak size and width depend on 
number, interactions, and proteolytic rate. This proteolytic size-control mechanism is consistent 
with the phenomenology of pseudo-pilus length control in the general secretion pathway of bacteria. 

PACS numbers: 87.16.A-, 87.16.dt, 64.75.Jk 



I. INTRODUCTION 

Living cells control the sizes of subcellular structures. 
Mechanisms of size control include molecular rulers for 
bacterial injectisome length [ij, measuring cups for the 
length of the bacterial flagellar hook counting for 
telomere ends |3[, and equilibrium energetics for actin 
bundle radius [Jl . Size control is also evident in the length 
of eukaryotic flagella [sj or the size of lipid rafts . 

Size control is challenging in bacteria because of the 
strong stochastic effects expected in such small cells. It 
is especially interesting how bacteria control the size of 
extracellular macromolecular assemblies, such as bacte- 
rial secretion systems and pili. In this paper, we investi- 
gate a novel length-control mechanism that may apply to 
the pseudopilus (here "ppilus") of the type-II secretion 
system (T2SS) of Gram-negative bacteria [3] . In the gen- 
eral secretory pathway, proteins are first secreted across 
the bacterial inner membrane by the Sec or Tat systems 
then across the outer membrane by the T2SS. The T2SS 
pushes folded proteins across the periplasm, and out a se- 
creton in the outer membrane, using an assembling and 
disassembling ppilus that is thought to function as a pis- 
ton or plunger Q. The ppilus assembles from the ener- 
gized inner membrane and spans the periplasmic space. 
The primary pilin subunit that assembles into the pseu- 
dopilus of the T2SS is variously called PulG in Klebsiella 
oxytoca, XcpT in Pseudomonas aeruginosa, or more gen- 
erally GspG (here "G") and is homologous to the PilA 
pilin of the type-IV pilus used in twitching motility in, 
e.g., P. aeruginosa or Myxococcus xanthus |8|. 

For a functioning T2SS the ppilus length should span 
the periplasm, which is approximately 21 nm across @ 
- or 85 G monomers in the ppilus structure [l^. The 
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ppilus is not normally seen outside the cell 0, hence 
the "pseudo" prefix, indicating an effective length-control 
mechanism. Overexpression of G leads to visible extracel- 
lular ppili (lll - [l3l | which rules out the (fixed-size) molecu- 
lar ruler or measuring cup mechanisms of size regulation. 

Underexpression of a minor pilin (XcpX [ll|, PulK 
[l3| . or GspK — here "K") leads to long extracellular 
ppili. So, it has been proposed that simple stoichiometric 
control applies to T2SS ppilus length control [13]. How- 
ever, stoichiometric control faces the inherent challenges 
of precisely controlled protein expression [l5l | . Further- 
more, while a single "stoichiometric" ppilus with a fixed 
size pool of G monomers can have a narrow length dis- 
tribution peaked at the pool size (TB], this is not true of 
multiple ppih sharing a common pool of G monomers. 
With multiple ppili stoichiometric length distributions 
are exponential, as seen in stochastic simulations with 
more than one ppih (see Appendix A). Since five to ten 
ppili are present on each individual bacteria [l7| , stoichio- 
metric length control of each individual ppilus requires 
an additional mechanism to partition G proteins equally 
between the ppili. 



II. MODEL 

The observations of G-G interactions [3l and of G- 
clusters in individual bacteria [l^ [l^ are consistent with 
clusters of G that could be associated with each ppilus, 
as illustrated by the dashed regions in Fig. 1. The size 
of each cluster, i.e. the number of G monomers, would 
then determine the maximal length of the associated 
ppilus, converting a ppilus length-control problem into a 
G-cluster size-control problem. Nevertheless, thermally 
driven evaporation, condensation, and diffusion (Fig. 1 
(e), (c), and (D) respectively) will destabilize sponta- 
neous partitioning of G among many clusters. The subse- 
quent coarsening of the size-distribution, treated by Lif- 
shitz, Slyozov and Wagner (LSW) [l^, would lead to a 
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FIG. 1: (Color online) Cartoon of the bacterial periplasmic 
space showing the inner (IM) and outer (OM) membrane, 
with pilin "G" proteins (red circles). Two pseudo-pili are 
shown, extending from their inner membrane base towards 
their outer membrane secretons (green squares) . The IM clus- 
ters associated with the ppili are indicated by the dashed 
circles. The number of G proteins in the ppilus together 
with its associated cluster is unchanged by ppilus assembly 
or disassembly (not shown). Numbers change by the physical 
processes of condensation (c), evaporation (e), proteolysis or 
recycling (-), and insertion (-I-). Diffusion (D) of monomers 
couples the clusters. 



single large cluster (condensed phase) of G in equilibrium 
with small clusters of G (vapour) in the bacterial inner 
membrane. The net growth of large diffusively-coupled 
2d clusters can be expressed in terms of their radius R 
by dR/dt = A (1/Rc — /R, where A is a constant. 
Cluster growth is the net result of a condensation flux 
(due to the supersaturation of G monomers in the inner 
membrane) and an evaporative flux (due to the curved 
boundary of the cluster within the membrane). Rc is the 
critical cluster size above which clusters grow and be- 
low which they shrink, and it typically grows with time 
due to the decreasing supersaturation associated with in- 
creasing average cluster radius via the Gibbs- Thompson 
effect. Growth of the average cluster size is also asso- 
ciated with a decreasing number of clusters, and with a 
broad distribution of cluster sizes [2^. This is qualita- 
tively unchanged for collections of smaller or irregularly- 
shaped clusters. 

In a biological system we should add both protein syn- 
thesis and proteolysis. Protein synthesis will simply con- 
tribute to the supersaturation, but proteolysis will add a 
new degradative term to the LSW dynamics proportional 
to the number of monomers in the cluster. Bacterial pro- 
teolytic mechanisms include cytoplasmic ATP-dependent 
proteosomes [2l|, ubiquitin-likc targeting systems [2^ . 
and periplasmic ATP-independcnt proteases [1^. Non- 
dcgradative recycling of components away from the mem- 
brane has an equivalent effect j24l - |26j . In terms of the 
number of monomers in a large 2d cluster of size N we 
will then have 



dN E 

—— = S = - aN, 

dt Jn 



(1) 



where S corresponds to condensation due to supersatu- 
ration, E corresponds to evaporation, and a is the pro- 
teolytic rate. Eqn. [T]is an approximate mean-field equa- 
tion, without stochastic effects or higher-order curvature 



corrections expected for smaller clusters. However, it il- 
lustrates (see inset of Fig. how intermediate values of 
the proteolytic term should generically stabilize coarsen- 
ing to a steady-state cluster size (red dot in inset), while 
for large enough proteolytic rate no stable fixed-point ex- 
ists (blue curve in inset). 

To test these ideas we model G clustering in the bacte- 
rial membrane with a stochastic Ising lattice-gas in two- 
dimensions (see e.g. [13] )■ The dynamics are (conserved) 
particle-exchange, subject to a Metropolis acceptance cri- 
terion with a reduced interaction energy J = J/{kBT). 
We supplement these dynamics with a dimensionless pro- 
teolysis rate a (per monomer per timestep), where (to 
minimize finite-size effects) monomers are removed and 
replaced at random positions in the system. Using a 
typical monomer size of Aa; ~ 5nm and diffusivit y o f 
D ^ 70000nm2/s = Sx'^/{4St) we have At « IQ-'^s [23, 



and the proteolytic lifetime is t = At /a. In units of 
Ax, we use a linear lattice size L ^ 400 and check that 
finite-size effects are not significant. A typical bacterium 
is larger in every direction, with L w 2000 [l^] . 



III. RESULTS 

Figure [H illustrates the steady-state results of the 
stochastic simulations. For a high proteolytic rate (a = 
10~^, or protein lifetime t « 10s) no large domains are 
seen. At an intermediate proteolytic rate (a = 5 x 10~*, 
or r w 2000s), stable domains with a characteristic 
size are seen. At a low proteolytic rate (a — 10~^, or 
T « lO^s) proteolysis is insufficient to stabilize the largest 
domain against coarsening to the limits of the system 
size. For the remainder of the paper, and for interme- 
diate proteolytic rates, we quantify the peak size Nq as 
well as the full-width-half-maximum W of the steady- 
state distribution of domain sizes, and explore how they 
vary with proteolytic rate, total membrane density of G, 
and reduced interaction J. 

The growth dynamics of a given cluster of size TV is 
given by the difference of incoming and outgoing flux, 
N = N+ (N) - iV„ (N) . This gives us the transition prob- 
abilities of monomer addition r+(A^) = N+/{N_ + N+) 
and subtraction r_(iV) = N_/{N_+N+). If P(7V) is the 
resulting steady state probability distribution of clusters 
of size N, then the detailed balance condition of the tran- 
sition probabilities is r+(iV- l)P(iV- 1) = r_{N)P{N). 
Approximating P{N) as a continuous distribution of 
cluster sizes (where dP/dN ~ P{N) - P{N - 1)), then 
in steady state we have 



P{N) ~ P(l)cxp 



N 



r+(n-i) 



dn 



(2) 



where we choose P{1) such that P{m) = 1. 

We have measured the steady-state evaporation and 
condensation rates N±{N), and used the resulting tran- 
sition probabilities T±{N) in Eqn. [2] to compare with the 
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FIG. 2: (Color online) Cluster size probability distribution, 
P{N), vs. number of monomers in a cluster, A^, for low 
{a = 10~^, green, right), medium (a = 5 x 10~*, red, cen- 
ter), and high (a = 10~^, blue, left) proteolytic rates (per 
protein per timestep, where At ~ 10~*s). Other parameters 
are J = 1.67 and p = 2 x 10"''. The dotted lines represent 
the reconstruction of the probability distribution via Eqn. [2] 
Intermediate proteolytic rates stabilize a peaked steady-state 
cluster size distribution with multiple large clusters, as illus- 
trated by the corresponding snapshot of the system. The in- 
set is a stability diagram from Eqn. [TJ illustrating the effects 
of low (green, top curve, with only an unstable fixed point), 
medium (red, center curve, with an additional stable fixed 
point as indicated), and high (blue, bottom curve, with no 
fixed points) proteolytic rates on the domain size dynamics 
dN/dt vs. domain size N . 



measured P{N). We find reasonable agreement (see e.g. 
the dotted lines in Fig.[2]or the crosses in Fig. [3]), showing 
that the dispersity of domain sizes around the stable fixed 
point No (where N+{No) = N^{No)) is due to stochastie 
fluctuations driven by the finite fluxes. Intuitively, and 
as proven in Appendix B, the stochastic width W of the 
size-distribution increases either with increasing absolute 
flux at the flxed point (i.e. |iV±|) or with decreasing net 
(stabilizing) flux (|iV+ — A^_|) near the stable fixed point. 

Fig. [3ja) shows the peak position and the width of the 
steady-state cluster-size distribution vs. proteolytic rate 
a. Increasing a decreases the stable domain size until, 
for a > 4 X 10~^, the distribution becomes peaked at 
A''o = 0. This is consistent with Eqn. [T] and the sta- 
bility plots in the inset of Fig. [2l The width decreases 
with a as the net stabilizing flux near the stable point 
increases. As shown by the inset of Fig.[3ja), the relative 
width W/Nq nevertheless increases with a. The fraction 
/ of G-monomers that are found in larger clusters (under 
the second peak in Fig. 2) slightly decreases with a but 



remains a considerable fraction (more than half) of the 
total, which indicates that proteolysis can be an efficient 
size-control mechanism. Experimentally, it appears that 
at least 20% of G-monomers are in ppili [1^ . 

In contrast, the steady-state peak cluster size Nq, the 
width W near the peak, and the fraction of G in the peak 
clusters stay roughly constant over a large range of aver- 
age G densities, p = Ng/L^, as shown in Fig. [3ljb). We 
understand this as an effective (nonequilibrium) coexis- 
tence between a fixed density of G monomers and excess 
G in clusters — increasing p simply increases the num- 
ber of clusters without significantly changing their size- 
distribution. (When approximately one cluster is seen 
in the system, for smaller p, finite-size effects do appear 
— decreasing both A'o and W . This is beginning to be 
apparent at the lowest p in Fig. [2Ib). Strong finite-size 
effects are seen in the condensed-cluster fraction /, but 
simply arise from an approximately constant vapour den- 
sity as p varies.) The lack of strong dependence of the 
cluster size A'o on the total membrane density p is ad- 
vantageous in terms of robust control of cluster size in 
the face of stochastic protein expression. 

At a fixed proteolytic rate and expression level, the 
effects of varying J/^ksT) (i.e. G-G interactions) are 
shown in Fig. [3Jc). Thermal evaporation decreases with 
increasing J — leading both to an increasing / (see inset 
of Fig. [3ljc)), a decreasing effective supersaturation, and 
a smaller stable Nq. The fractional width W/Nq of the 
peak of the cluster size distribution is narrower for weaker 
interactions, but at the same time a smaller fraction of 
G-monomers are in clusters. 



IV. DISCUSSION 

We find that intermediate levels of proteolysis controls 
the natural coarsening instability of condensed clusters, 
and leads to steady-state clusters with a relatively nar- 
row distribution of sizes. For an intermediate proteolytic 
rate a = 5 x 10^* (turnover time r ~ 2000s), we obtain 
a cluster size iVg ~ 80 — remarkably close to the 85 G 
required to assemble a ppilus that spans the periplasmic 
space [1, The fractional width is then about 30%, 
consistent with the lack of extracellular ppili observed 
under normal conditions. For cluster size-control to ef- 
fectively control ppilus length, we predict that G-clusters 
are associated with the secreton base of the T2SS. Mu- 
tational variation of G-monomers to affect either their 
proteolytic susceptibility {a, see e.g. Fig. [2Ia)) or G- 
G interactions (J, see Fig. ^c)) should also affect the 
distribution of ppilus lengths. We believe stoichiometric 
mechanisms, via GspK |lll - [l3l |. control the ppilus length 
from individual associated G-clusters, while our prote- 
olytic size control mechanism ensures that multiple G- 
clusters remain approximately equally sized. 

We have shown that a novel mean-field fixed-point in 
the cluster size distribution arises from proteolysis, while 
the dispersity of cluster sizes around the fixed point arises 
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FIG. 3: (Color online) (a) Dependence of the steady-state most-likely cluster size A^o (red circles) and the full-width-half- 
maximum W (green squares) of the cluster-size distribution vs. the proteolytic rate a. Crosses (x) indicate corresponding 
values derived from fluxes near A'^o "via Eqn. (2] The inset shows the relative width W/No (blue squares) and the fraction of 
G-proteins in large domains (purple circles) vs. a. Error bars indicate statistical errors, which are often smaller than point 
sizes. Other parameters are J = 1.67 and p = 2 X 10~^ (b) The top fi gure shows size A^o (red circles) and width W (green 
squares) vs. average membrane density p = Ng/L^ . The bottom figure shows the dependence of relative width W/Nq (blue 
squares) and the cluster fraction / (purple circles). Other parameters: J = 1.67 and a = 5 x 10~*. (c) Size A^o (red circles) and 
width W (green squares) vs the clustering interaction J = J/iksT). The inset shows the relative width W/No (blue squares) 
and the fraction of G-proteins in large domains / (purple circles) vs. J. Other parameters: p = 2 x 10~^ and a = 5 x 10~*. 



from stochastic gro'wth and shrinkage of clusters. The 
noise associated with proteolysis is intrinsically multi- 
plicative, in that proteolysis only targets existing G pro- 
teins. Our latticc-gas model naturally implements both 
proteolysis and the thermal evaporation and condensa- 
tion of clusters in a membrane. 

The result is a monodisperse cluster-size distribution 
with a non-zero peak size, which qualitatively differs from 
some earlier work on lipid raft sizes in membranes that 
only found distributions peaked at iVg « [1^ [26| . 
We believe this is due to the approximate evapora- 
tion/condensation dynamics [l^ or the additive noise 
[26| used in those works. In contrast, earlier coarse- 
grained models of ternary mixtures with recycling [2^ 1 — 
also applied to lipid nanodomains — did recover a non- 
zero peak size, though did not include recycling noise. 
Our microscopic two-component model is simpler, and 
the mean-field flow we present shows how the results are 
expected to be generic for proteolysis or recycling. 

Proteolysis is not just for cellular cleanup. Targeted 
degradation can adjust timescales and levels of transcrip- 
tion or translation (see, e.g. [lH). We have shown how 
it can also be used to control cluster sizes within the cell. 
Proteolysis contributes a new "evaporative" term that 
limits coarsening with a novel mean-field fixed point for 
the cluster size (inset of Fig. [5]). The same mechanism 
will qualitatively apply whether the proteolysis targets all 
proteins in a cluster (as in Eqn. [1]), appropriate for cyto- 
plasmic or periplasmic proteases — or whether it targets 
the cluster periphery as might be appropriate for mem- 
brane associated proteases. In both cases the loss term in 



Eqn. [T] will grow with A'', and so will lead to a stable fixed 
point at some iVp. Proteolysis provides a new size-control 
mechanism to cells. We expect that proteolysis or anal- 
ogous degradation terms, such as recycling, are widely 
exploited to achieve monodisperse steady-state clusters 
in other biological systems. 
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Appendix A - Stoichiometric length-control 

We present here analytical limits as well as computa- 
tional simulations for a stoichiometric model for length- 
control. This system consists of M polymerized pili shar- 
ing a common pool of Ng = MNp monomers, where Np 
is the average number of monomers per pilus. Each pilus 
has a probability of growing and p_ of terminating 
growth and completely disassembling, and these prob- 
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FIG. 4: (Color online) Probability distribution Pm{1) of the 
maximum length for various number of pili (M = 1, 2, 3, 4, 6) 
sharing the same common pool of monomers. We use t = 0.3 
and Np — 40, though qualitatively similar results are seen 
with other parameters. The continuous lines are the result 
of the stochastic simulation. The red dashed lines represents 
the two exact limiting distributions pi and poo, overlaying the 
M — 1 data and close to the M = 6 data, respectively. 



exponential distribution 

P^H) = (4) 

For the intermediate regime, with finite Af > 1, we per- 
formed stochastic computer simulations — as illustrated 
in Fig. m The exponential limit is quickly approached 
for M > 4. However, even for M = 2 the peak around 
A^r, seen for A/ = 1 is lost. 



Appendix B - Analytical monotonicities of the width 

For the proteolytic size-control mechanism described 
in the text, we detail here the formal derivations of the 
variation of the full width half maximum (FWHM) W of 
the non-zero peak in the size distribution. We consider 
the evaporative and condensation fluxes J+ and J_, re- 
spectively, close to the stable fixed point size N = Nq, 
where J+ = J- = Jq (see Fig. [5]). We allow for linear 
dependence of the fluxes near the fixed point 

J+ = a{N-No) + Jo (5) 
J_ = b{N-No) + Jo (6) 

where the stability of the fixed point requires J_ > ,/+ 
for N > Nq and J+ > J_ for N < Nq, corresponding to 
the requirement that b < a. 



abilities are dependent only on the pool of monomers. 
The total number of monomers is constant, so that Ng = 
n + '^f^^ li, where li are the lengths of the pili and n the 
number of monomers remaining in the monomeric pool. 
Since the growth rate is proportional to n, while the ter- 
mination rate t is fixed, then the growth and termination 
probabilities are p+ (n) = n/{n+t) andp_(n) = t/{n+t), 
respectively, where + ]3_ = 1 . 

For a single pilus, the probability of achieving a maxi- 
mum length / can be written as the product of the indi- 
vidual probabilities: pi{l) = p+(A^g) x p+{Ng — 1) x • • • x 
p+{Ng 1)) X p_(iVg - I). This leads to 



Pi{l) 



tT{Ng + l)T{Ng +t-l) 
T{Ng+t+l)T{Ng-l+l) 



(3) 



which is peaked around Np for t < 1, as shown with M = 
1 in Figure |4l A single pilus can achieve length-control 
by assembling most of the available pool of monomers 
before disassembly is triggered. 

For the more biologically appropriate case of multiple 
pili per cell, we can use a mean-field approximation. In 
this case, t becomes Alt, and Ng becomes MNp — Ml 
where / is the average length of the pili, and thcnpA/(0 ~ 
p+{MNp - Ml) X p+{MNp -Ml-l)x---x p+{MNp - 
Ml - {I - 1)) X p^{MNp ~ Ml- I). The p+ are constant 
at 0(1/Af), so that in the limit M ^ 1 we recover an 
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FIG. 5: Figure illustrating the mean-field fixed-point cluster 
size A'o that occurs where the evaporative and condensative 
fluxes are equal (J_ = J+ = Jo). The fixed-point is stable 
when evaporation is stronger than condensation (J_ > J+) 
for N > A'o — as illustrated. The demonstration of mono- 
tonicities in this appendix are for N ~ Nq, where a linear 
approximation for J± holds. 



These fiuxes give us the transition rates 

an + Jo 



cn - 
bn 



2 Jo 

-Jo 



2 Jo' 



(7) 
(8) 
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where we define n = N — Nq and c = a + b. The detailed 
balance condition is p{n) = p{n — l)r_|_(7i — l)/r_(n). 
We define 



A{n) ^ n 



m— 1 

cn/2 



r+(m-i) 
r_(m) 



^0 mil ('"a 



(9) 



6n 



•^0 nriiM+^o) 



A monotonically decreases with n, since r+(n — 
l)/r_(n) = (a(n - 1) + Jo)/(6(n - 1) + Jq) < 1 where 
h> a. The FWHM condition, A{W/2) = 1/2, allows us 
to determine how W must respond to changes in Jq and 
A EE &- a. 

Varying Jg at FWHM we have 



where c/2 < h and a < 6 so that dlogA/dJo > and 
A monotonically increases with Jq. We conclude that W 
increases with increasing Jo- 
in the same manner. 



dlog^ 



W 



W 



WA + 2aW + 4 Jo T^A ^ aW + 2 Jo 



nA + no + Jn 



(11) 



dlogA 1 
dJo "cVF/4 + Jo 

W/2-1 

+ E 

n=l 



5W^/2 + Jo 
1 1 



na + Jo n& + Jo 



(10) 



is negative, so that A monotonically decreases with A = 
h — a. We conclude that W decreases with increasing 
6 — a. 
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